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We discuss the problem of adding random matrices, which enable us to study 
Hamiltonians consisting of a deterministic term plus a random term. Using 
a diagrammatic approach and introducing the concept of "gluon connected- 
ness," we calculate the density of energy levels for a wide class of probability 
distributions governing the random term, thus generalizing a result obtained 
recently by Brezin, Hikami, and Zee. The method used here may be applied 
to a broad class of problems involving random matrices. 
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1 Introduction 

Some four decades ago, Wigner [1] proposed studying the distribution of 
energy levels of a random Hamiltonian given by 

H = ^ (1) 

where cp is an N hy N hermitean matrix taken from the distribution 

Zj 

with Z fixed by / d(pP{ip) = 1. This problem has been studied intensively by 
Dyson, Mehta, and others over the years [2,3,4]. Two years ago, Brezin and 
Zee discovered that, remarkably, while the density of eigenvalues depends 
[5] on V, the correlation between the density of eigenvalues, when suitably 
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scaled, is independent [6] of V. This universality was also obtained earlier 
and independently by Ambj0rn and collaborators [7]. Since then, it has been 
clarified and extended by other authors [8,9,10], studied numerically [11], and 
furthermore, shown to hold even when the distribution (Q) is generalized to 
a much broader class of distributions [12]. We expect that the discussion to 
be given below will hold also for this broader class of distributions, but for 
the sake of simplicity we will not work this through here. 

In recent work [13], Brezin and Zee have generalized this Wigner problem 
to the case of a Hamiltonian given by the sum of a deterministic term and a 
random term 

H = Ho + ip (3) 

Here Hq is a diagonal matrix with diagonal elements e^, i = 1,2,...N, and 
ip a random matrix taken from the ensemble (^. For the Gaussian case, 
namely with V{(p) = ^93^, Pastur [14] has long ago determined the density of 
eigenvalues. The work described in [13] went beyond Pastur's work in that 
the correlation function between the density of eigenvalues in the Gaussian 
case was also determined. More recently, in a work with Brezin and Hikami 
[15], we managed to determine the density of eigenvalues for V{ip) = + 
gif"^ to all orders in g. The correlation function was also computed, but only 
to first order in g. 

This problem of "determinism plus chance" may be regarded as a generic 
problem in physics, and as such represents a significant generalization of 
Wigner's problem. For example, consider an electron moving in a magnetic 
field and scattering off impurities. We note that these "deterministic plus 
random" problems are considerably more difficult than the purely random 
problems defined in (|l|) and (H). A standard approach to solving the purely 
random problem involves diagonalizing the random matrix (p and then use 
orthogonal polynomials to disentangle the resulting expression. Clearly, in 

we cannot diagonalize if without un-diagonalizing Hq and thus the or- 
thogonal polynomial approach fails. 

In this paper, we point out that the problem given in is a special case 
of a broader class of problems involving the addition of random matrices. 
The deterministic Hamiltonian Ho may in turn be replaced by a random 
Hamiltonian. Indeed, a deterministic matrix is but a special case of a ran- 
dom matrix. We will extend the work of Brezin, Hikami, and Zee [15] and 
determine the density of eigenvalues for the Hamiltonian given in (^) for an 
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arbitrary V. 

Our work is inspired by recent advances in the mathematical hterature 
involving the theory of non-commutative probability and operator algebra 
[16,17,18]. A number of physicists have already brought these advances to 
the attention of the physics community [19,20,21]. While our work is thus 
inspired, we will not be using the mathematical approach given in [16,17,18], 
but instead will be based on the diagrammatic approach developed in [13] 
and subsequent work [22,15]. 

2 Adding Random Matrices 

Consider a Hamiltonian given by 

H = (pi + (p2 (4) 
with the matrices ipi^2 taken from the probability distribution 

P{^i, ^2) = ie-^*^[^^(^^)+^^(^^)l = Pi(^i)P2(^2). (5) 

Notice that the probability distribution factorizes. (This is known as "free" 
in the mathematical literature.) The problem defined in represents a 
special case. Previously, with D'Anna and with Brezin we have studied the 
problem given in (1) but for the more difficult case [23,24] in which P{ipi, (p2) 
contains terms linking ipi and ip2- Indeed, a detailed determination of the 
correlation function over all "distance scales" is a non-trivial problem even 
for Gaussian distributions [24]. The discussion in this paper goes through 
precisely because (pi and ip2 do not couple to each other in P(v?i, ip2)- 

Let us now mention a few necessary definitions. Define the Green's func- 
tion 

The density of eigenvalues is then given by p{jj,) = (^-^tr6{jj, — H)'^ = — ^ImG(/x+ 
ie) . In this paper we focus on the density of eigenvalues, leaving the correla- 
tion for a future work. Note that the factors of are chosen in our definitions 
such that the interval over which p{fi) is non-zero is finite [i.e., of order N^) 
in the large A^ limit. 
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We may regard the distribution as defining a (0 + 0)-dimensional field 
theory. The Feynman diagram expansion is then simply obtained by expand- 
ing G{z) in inverse powers of z and doing the integrals in (^. As explained 
in [13], it is useful to borrow the terminology of large N quantum chromo- 
dynamics [25] from the particle physics hterature, and speak of quark and 
gluon lines. See figure (1) for a graphical representation. (It is of course not 
necessary to use this language, and readers not familiar with this language 
can simply think of the diagrams as representing the different terms one en- 
counters in doing the integral in (|])). The quark propagator simply comes 
from the explicit factor of z in and is represented by a single line and given 
by p The quadratic terms in Vi(v?i) and V2{(p2) determine the gluon propa- 
gator, represented by double lines. Here, in a minor departure from large N 
quantum chromodynamics we have two types of gluons, corresponding to (pi 
and ip2. The gluon propagators are proportional to 

(viAi) oc So^pSis';^. (7) 

The non-Gaussian terms in Vi{ipi) and V2{(p2) describe the interaction be- 
tween the gluons. 

The reason that we can solve this problem is because, while the two types 
of gluons have arbitrarily complicated interactions among themselves, they 
do not interact with each other. Note that while the gluons both interact 
with the quark, our problem is such that we do not have to include quark 
loops and thus the quark does not induce interaction between the two gluons. 
This is clear from the definition of our problem. Another way of saying this 
is to note that the Green's function may be represented, by using the replica 
trick, as 

G{z) = Urrin^o J Dt/j^ Di/jDipP{ip)ijliJie- S^=i ^-(^-'^)'/'- (8) 

Note that in this language the ip^s represent the quark fields and cp the gluon 
fields. The interaction between gluon and quarks are given by ipafi'a- (Color 
indices are suppressed here.) The interaction of the gluons with each other 
is determined by P{(p)- Since internal quark loops are proportional to the 
number of replicas n, they vanish in the n ^ limit. 

Let us then calculate the Green's function, which as usual can be written 
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as (see figure 2) 



G{z) 



z-j:(z) 



(9) 



in terms of tlie one-particle irreducible self energy = The 

self-energy is then determined by the set of diagrams in figure (3) with the 
corresponding equation 



< j^trifi + < J^t^fl >ac G{z)+ < ^trifl >g^ G{zf + .... + (1^2) 

E < ^tr^t >.c Giz)'-' + (1 - 2) 
fc=i 



1 
G 
1 

G 



1 

< — tr— 
N 1 



+ (1 



+ (1^2) 
^ 2). 



In order to write this equation, we have to invoke the factorization of P{^pi, ^P2) 
which tells us that the two kinds of gluons do not interact, and the large N 
limit, which tells us that the two kinds of gluon lines cannot cross. 

We are led to introduce in (|1^) the notion of "gluon connectedness," 
denoted by "gc" henceforth. The necessity for this notion is illustrated by 
the shaded blob describing the interaction of the gluons in figure (3d): it 
should not include the diagram shown in figure (4): this class of diagrams is 
already included in figure (3b). In other words, a gluon connected blob with 
k external gluon lines is such that it cannot be separated into two blobs, with 
ki gluon lines and /c2 gluon lines respectively, (with ki + k2 = k of course). 
In the last line we have defined the "gluon connected Green's function" 



(10) 



Ggci{z) =< -^Jr 



N z-ipi 



> 



(11) 



and similarly Ggc2{z). The operations implied in (|IT]) are clearly allowed 
since < j^tr- >gc is a linear operation. Note also that we have not assumed 

In particular, we include a 
'gc in (0). 



that Va is an even function of its argument, 
possible tadpole term indicated by < jjtnpa > 



We should emphasize that the shaded blobs include interactions between 
gluons to all orders. It is very complicated, if not hopeless, to calculate these 
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blobs in terms of Vi and V2, but fortunately, as we will show below, we do 
not have to calculate them exphcitly. In our previous papers, we regarded 
the cubic, quartic, and so on, terms in Va as interactions and proceeded to 
calculate the Green's function and correlation function in terms of the various 
coupling constants. We follow a different strategy here, and try to express 
the Green's function G{z) directly in terms of Gi{z) and 6*2 (-2) where 

for a = 1,2 are the Green's functions for two separately and purely random 
problems. (The average in ( ]T2| ) is performed with the distribution Paifa) = 
_i_^-NtrVa(iPa) Qf course.) In this way, we attempt to bypass having to deal 
with Vi and V2 altogether. 

To see how to do this, let us go back to the simpler problem defined by (|^) 
and (|^). Following the same diagrammatic analysis leading to (|TUp we find 
that the Green's function G{z) and self energy for this simpler problem 
are related by 



(13) 



with, evidently. 



Goc{z) ^tr-l- (14) 



N z-ip 
Combining ([1^) and (^ we find 

^GU^) = z (15) 

For the sake of convenience, we may, with due respect to Green, somewhat 
fancifully define a "Blue's function" by 

B{z) ^ ^G,,(^) (16) 



Thus, we learn from (|l^) that the Blue's function is the functional inverse of 
the Green's function 

BiGiz)) = z (17) 

From the normalization of the probability distribution P{f) we obtain triv- 
ially the "sum rule" G{z) ^ - as z ^ 00, thus implying that the Blue's 
function B{z) ^ - as z 0. 



6 



Let us now go back to the more involved problem defined by (^). First, 
we define for a = 1,2 two Blue's functions Ba as the functional inverse of 
Ga respectively. We now see that (PDD, when combined with says simply 
that 

z + ^ = B,{G) + B^iG) (18) 
Thus, the law of addition for the Blue's function is given by 

Bi+2iz) = Biiz) + B2iz) ~ ^ (19) 

This equation tells us how to obtain the Blue's function associated with 
fi + V^2 from the Blue's functions associated with ipi and ip2- 

The procedure for determining the Green's function and hence the den- 
sity of eigenvalues of the problem defined by (|l|) and @ is then as follows: 
given the Green's functions Gi and G2, determine the corresponding Blue's 
functions Bi and B2 by functionally inverting Gi and G2 respectively, cal- 
culate -Bi+2 according to (pJ]), then determine the functional inverse of -B1+2 
to find the desired Green's function G{z). 

Let us remark briefly on the connection to the mathematical literature. 
Voiculescu [16] has introduced the "/^-transform." It turns out that the R 
function discussed by mathematicians is simply related to B hj B{z) = ^ + 
R{z). In fact, we see that the Dyson-Schwinger equation (|^) when combined 
with ( p!?!) gives simple B{G{z)) = + S(z). Thus, the R function of the 
mathematicians is nothing but the self-energy E of the physicists expressed 
in terms of different arguments: 

R{G{z)) = S(z) (20) 



3 Addition Algorithm at Work 

Let us now proceed by building up from a few simple examples. In the most 
trivial case, (p is not random at all, but fixed to be a constant c times the 
unit matrix. Then from the Green's function G{z) = we find the Blue's 
function B{z) = c + K For a slightly less trivial example, let be a diagonal 
matrix with matrix elements given by Cj with i = 1,...,N. The Green's 
function is given by 
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Then the corresponding Blue's function is determined by 



Next, let P{f) be Gaussian (that is, V{ip) = tr^ip"^). Then as is well known 
(see for example [13]), the Green's function is determined by 

In other words, 

1 



G{z) = -{z - ^/^^) (24) 
Substituting z ^ B into (0), we obtain immediately that [26] 

B{z) =z + ^ (25) 

Thus, in this simple case, the Blue's function B, which of course contains 
the same information as the Green's function G, actually has a simpler form 
than G. 

Now we are ready to do our first non-trivial problem. Consider the prob- 
lem defined in (^. Since we know from (p^) and (25) the Blue's functions 
corresponding to the two terms in the Hamiltonian, we learn immediately 
from (^) the Blue's function for H: 

Bi+2{z) = Bi{z) + z + ^ - ^ = Bi{z) + z (26) 



with Bi determined by ( P^ with the substitution B Bi. The desired 
Green's function G{z) is now determined by solving for the functional inverse 
of the function -B1+2, that is, by the equation 

Si+2(G) = z (27) 

or equivalently, upon using (p6l), 

Bi{G)=z-G (28) 
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Let us now evaluate the two sides of this equation with the function G'i(-). 
Since Gi{Bi{G{z))) = G{z) we obtain immediately 

G{z) = Gi{z - G{z)) (29) 

Noting that Gi is given by (|2l|) with the substitution G — > Gi, we have 

precisely the classic result of Pastur which was obtained diagrammatically in 
[13]. 

After these simple exercises, we can now immediately go on and solve 
the general version of the problem defined in @: find the density of energy 
levels of a Hamiltonian given hj H = Hq + (f with (f drawn from the general 
distribution @. With a slight shift in notation, let us call the Blue's function 
associated with Hq and with ip respectively Bq and B2. Then the Blue's 
function associated with H is given hjB{z) = Bo{z)+B2{z)--. Substituting 
in this equation z G{z) (where G{z) is the unknown Green's function 
associated with H), we find immediately that 

Bo{G) = z + ^ - B2iG) (31) 

Anticipating the next step, we define 

S(z) = B^iGiz)) - ^ (32) 



With this definition, we write ( [3T| ) as 

Bo{G) = z -J:{z) (33) 

Let us now evaluate both sides of (|31| ) with the Green's function Gq{-) asso- 
ciated with Hq. We find instantly that 

Go {Bo{G{z))) = G{z) = Go{z- S(z)) = 1 ^ z-e,^-E{z) ^^^^ 

Let us repeat this trick: rewrite ( ^2]) as B2{G{z)) = T,{z) + and evaluate 
both sides with the function G2{-). We obtain 

G{z) = G2{m + ^). (35) 
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These two equations, ( |34D and (pSj), allow us to determine the two un- 
known functions G{z) and ^{z), provided we know the Green's function 
G2{z). But what is the Green's function G2{z)7 It is just the Green's func- 
tion associated with the random matrix ip drawn from the general distribution 
(H). But this was obtained by Brezin et al [5] almost twenty years ago. These 
authors told us that (for V{z) an even polynomial for the sake of notational 
simplicity) 

G^iz) = ^[V'iz) - Piz)V^^^] (36) 

Here V'{z) = P{z) is a polynomial, and a determines the endpoints of 
the spectrum of eigenvalues. The quantities P{z) and a are determined [27] 
by the "sum rule" that G2 —>■ - as z ^ oo. 

In summary, and repeating various equations for clarity, we have obtained 
the following result. For a Hamiltonian of the form H = H^ + ip with the ran- 
dom matrix Lp drawn from an arbitrary distribution defined by V{ip) (taken 
to be even for simplicity), we can determine the Green's function G{z) and 
hence the density of eigenvalues by solving simultaneously the two equations 



GW = G2(E(j) + — -), (38) 



and 



where the function G2 is given by (pGf). In general, for an arbitrary set of e^'s 
and a non-Gaussian V , these equations can only be solved numerically. 

It is clearly of some notational benefit to give the combination appearing 
in ( |38D a name: a{z) = T,{z) + We can then simplify ( pS]) shghtly to 
[28] 

P\a){a^-a^) = {V'{a)-2Gf (39) 

Thus, we can use (^) to determine a, and hence S, in terms of G. Plugging 
this into (|37| ) then gives us an equation for G. 

As mentioned earlier, Brezin, Hikami, and Zee [15] recently used the equa- 
tion of motion method and a detailed diagrammatic analysis to determine 
the Green's function for the problem in (^ with the distribution defined by 
V{ip) = + gip"^. It is straightforward, although slightly tedious, to verify 
that for this simple case, (^9|) reduces to equation (4.15) in [15]. The analysis 
given here is considerably simpler. 
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4 Energy Bands 



The function B contains the same amount of information as G, since one 
function is the inverse of the other. As an example, consider the problem of 
determining the endpoints of the energy spectrum. Near the endpoint, call 
it a, the density of states p(/i), that is, the imaginary part of G{z), vanishes 
generically like (a — /x)2. Thus, the endpoint is determined by the equation 

^l„. = oo (40) 

Now consider the defining equation for B: namely G{B{w)) = w. Differen- 
tiating, we obtain 

dG dB / Ai\ 

dBd^ ~ ^ ^ 

Thus, we obtain an alternative equation for determining the endpoints a of 
the energy spectrum of a random Hamiltonian: solve the equation 

dB, 
dw 

Let us now apply these simple considerations to determine the endpoints of 
the spectrum of a Hamiltonian of the form H = Hq + ip. 

As remarked above, with an arbitrary Hq (that is, an arbitrary set of e^'s) 
and a general V, we can hardly expect to solve (0) and ( p^ ) analytically. 
Let us retreat to the case in which (f is taken from a Gaussian distribution, 
so that in (0) can be replaced by G{z). (In other words, the solution 
of (|38D is simply S(z) = G{z).) Thus, we have to solve 

Gw^Ej^Trbw (43) 

We will be interested in the case in which the Hamiltonian Hq exhibits de- 
generacy. Let K be the number of distinct ej's (with K < N of course.) To 
determine G we have to solve a polynomial equation of degree K + 1. 

As explained above, in general, given an equation determining G{z), we 
simply substitute z — > B{w) and G{z) — > w into that equation to obtain the 
equation for determining B{w). Thus, in the present example, the function 
B is determined by 

1 ^ 1 
w = —y — — ^ (44) 
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We note that this is a polynomial equation of degree K, that is, of degree 
one lower than the equation for determining G. In particular, for the case 
of K = 2, while we have to solve a cubic equation to determine G, we only 
have to solve a quadratic equation to determine B\ 

We will now exploit this fact and study the Hamiltonian H = Hq + (p for 
= 2, in other words, with no loss of generality we take the symmetric case 
with the deterministic piece 



^^0= " (45) 




and take to be a Gaussian random perturbation. We can imagine some 
possible physical applications of this Hamiltonian. For instance, consider 
electron scattering on impurities in a spin dependent quantum Hall fluids. 
In the absence of impurity scattering we have a spin up Landau level sepa- 
rated in energy from the spin down Landau level, with a Zeeman splitting 
of 2e. Impurity scattering is represented schematically by A model not 
precisely of this type, but in the same spirit, was proposed by Hikami, Shi- 
rai, and Wegner [8] and has been studied by a number of authors [9,10]. The 
reader can no doubt concoct other possible situations represented, at least 
schematically, by a Hamiltonian of the type considered here. 

For e = 0, the density of states of H is given by the familiar semi-circle 
law p{^) = -^y/A — which we obtain easily from (3). Note that the half 
width of the spectrum is equal to 2. In the other limit, e >> 2, the density of 
states decompose into two pieces. At some critical value ec the two disjoint 
pieces in the density of states touch and merge into one piece. An interesting 
question is whether ec is larger or smaller than 2. A simple physical argument 
based on level repulsion would suggest that ec < 2. 

The equation (|43|) for G becomes 

While an explicit solution of this cubic equation may be written down, it 
is rather unwieldy. In contrast, as just explained, B satisfies the quadratic 
equation 

zB{zf - {2z^ + l)B{z) + + (1 - e^)z = (47) 
whose solution can of course be immediately written down. 
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We would now like to determine the endpoints of the spectrum of H using 
the algorithm we developed above. We employ the following simple trick. 
Define the operator ^| with the vertical bar indicating that the derivative 
should not act on B. This operator is clearly useful in light of (^21). We now 
act on (^7[) with the operator (-2^1 — !)• The net effect of this operator is to 
simply replace z"^ in ( ^Tf ) by (n — We thus obtain a linear equation for 
B: 

{2z^ - 1)B -2z^ = (48) 



We solve this equation and (or somewhat more simply, the equation 
obtained by acting with ^| on (|47|)) simultaneously for z and B. The values 
of B thus obtained are in fact the values of a that we are trying to determine. 
Proceeding in this way, we find easily that the four values of a are given 

by 

1 (4e^-l + v/8i^)f ^^g^ 



2 V2e VS^^TT - 1 

(50) 



1 (4e2 - 1 - 



and — a+ and — a_. 

As a first check, we take the limit e ^ 0. We find that as expected, two of 
the a's, namely ±a_ become complex and hence unphysical. Indeed, a+ — 2, 
and we recover Wigner's classic result for a Gaussian random Hamiltonian. 

In the opposite limit with e large, we find 

a/2 

a+ = e(l + ^ + ...) (51) 

and 

a_ = e(l--i- + ...) (52) 

Thus, in the density of states the width of each of the two disjoint pieces, or 
bands, is equal to \/2 + = Thus, even when the two bands are far 
apart, the width is reduced compared to the Gaussian bandwidth by 

bandwidth 3 ^^^^ 



Gaussian bandwidth A^/2 
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Finally, the phase transition when the two disjoint bands merge into one 
occurs when a_ vanishes and becomes complex. This occurs at 

ec = 1 (54) 

As expected, ec is less than the unperturbed Gaussian half- width of 2. Inter- 
estingly, level repulsion reduces the naive expectation by precisely a factor 
of two. 

Note that for the purpose of this discussion we never had to solve for B 
explicitly, but the solution of (0) is of course easy enough to write down for 
the records: 

2z^ + l + Vl+46^^2 
B{z) = — (55) 

The function B{z) has two branch cuts starting at Zc = We check 

easily that the addition law in ([TD|), giving B{z) = Bq{z) + z, with Bq the B 
function associated with Hq of course, produces the same result. 

Now that we have discussed a specific example, let us now go back to the 
general case and write (36) as 

i,^^ (56) 

B — fi — z 

where (T(/i) = j^J^i^ifJ' ~ q) is the density of states of the deterministic 
Hamiltonian Ho. Again, with no loss of generality, we have set = 1 in the 
Gaussian distribution governing (p. The endpoint a of the spectrum of H is 
given by the value of B that simultaneously solves (|56D and the equation 



Note that as explained before (|57| ) is obtained by acting with the operator 

i\ on m- 

Given any cr(yu) we can thus determine in principle the endpoints of the 
energy spectrum of H. We can for example consider the limit of large ran- 
domess, that is, when the spectrum of Hq is small compared to the scale of 
the randomness. Expanding the denominators in (^) and (pTj) and defining 
< e^*^ >= / dfia{^)fi'^^ (assuming for simplicity that 0"(/i) is even) we obtain 

1 <e2> <e*> 
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and 

1 3 < > 5 < > 

^=(b37) + (1^+(1J3T)«+ (59) 

In the limit Hq or large randomness, we have z = = ±1 (and thus 
the endpoints a = ±2, as expected.) It is easy to solve these equations to 
any desired power. We obtain to 0{Hq) that the width of the spectrum is 
given by 

a = 2+ < > < >2 + < > +.... (60) 
A simple check shows that indeed, if we expand a+ in (|49| ) we find 

a+ = 2 + e^-^e^ + (61) 



This is consistent with (^) since in this simple model < >^=< >= e 
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5 Analytic Structure 

For the problem studied here B is simpler than G, but in other problems B 
is unfortunately more complicated. For instance, for a much studied class of 
probability distribution, the so-called trace class, defined by (^, and with 
V taken to be an even polynomial for simplicity, G{z) is given by [5], as 
mentioned earlier, 

G{z) = ^ [V'iz) - P{z)Vz^-a^] (62) 

where the polynomial P and the endpoint a are determined completely by 
the condition G{z) — > 1/z as z — > oo. Applying our rule of substituting 
G ^ w and z B, we obtain the equation 



2w = V'{B) - P{B)^B^-o? (63) 

which determines B. 

For V of degree 2s it is easy to see that B satisfies a polynomial equation 
of degree 2s — 1. Indeed, moving V'{B) in (|63D to the left hand side and 
squaring, we can rewrite (^3]) as 

wV'{B) -V? - Q{B) = (64) 
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where 4:Q{B) = V'{B)'^ - P{BY{B'^ - a?) is a polynomial of degree 2s - 1 
in B. (This is most easily seen by noting that G{B) 1/B for large B by 
definition.) 

We would like to study the analytic structure of B but we have not been 
able to make general and complete statements. One limited statement applies 
to the trace class just discussed. It is easy to show then that the branch cuts 
of B{w) is of the square root type. Let B{w) has a cut starting at w* and 
write B* = B{w*). At the tip of the cut, 

dBiw) , 

— ; U* = (65) 

aw 

or 

dG{B) , 

-J^\b.=^ (66) 

Note that these equations are the duals of the ones in (|40D and (|4^). Ex- 
panding, we find immediately that 



B = B*+\-^\bA A{w-w*)^ + .... (67) 




For the case studied in the previous section with Hq given as in (|45D we 
see from the explicit form given in (^) that B{w) indeed has a square root 
cut starting at w* = in agreement with the general analysis given here. 



6 Non-abelian Central Limit Theorem 

Gauss proved that if we add K random numbers Xj, i = 1,2,...K, with 
Xi taken from the probability distribution Pi{xi), then the normalized sum 
s = follows the Gaussian distribution in the limit K tending to 

infinity. This result plays an important role in physics and mathematics and 
accounts for the ubiquitous appearance of the Gaussian distribution. 

What if the variables Xi do not commute? In particular, suppose that 
instead of real numbers Xj, we have by random matrices <y9j, i = 1,2, ...K, 
taken from from the probability distributions 

P,((^.) = le-^*'-^'('^'). (68) 
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Does the normalized sum of all these matrices ^pt = J2i follow the 
Gaussian distribution in the large K limit? Intuitively, it seems that this 
ought to be the case. (Note that there exists another commonly considered 
class of random matrices, in which the element of the random matrices is 
each taken from a probability distribution (the same for each element). For 
this class, which we refer to as the Wigner class in our earlier work [12], the 
proposed theorem follows immediately from the usual abelian central limit 
theorem. Here we are speaking of the trace classes defined by ([681) . 

As it turns out, it is not difficult to generalize one of the standard proofs 
of Gauss's theorem to matrices. We will give this proof below. However, it 
would seem that the algorithm developed here for adding random matrices is 
almost tailor made to address this question of whether ipx = J2i follows 
the Gaussian distribution. It is mildly amusing to see how the Wigner semi- 
circle law emerges naturally. 

As we mentioned earlier, when we add two random matrices, in general 
it is difficult to determine explicitly the resulting G{z) for the sum of the 
two matrices. What we hope for here is that the large K limit will bring 
considerable simplification. This is indeed the case. To keep the formulas 
simple, let us again take Vi to be even. We expect that our conclusions can 
be easily generalized. In this case, Gi{z) is an odd function and hence Bi{z) 
is also an odd function. Thus, we write Bi{z) = ^ + zbi{z) with bi{z) an even 
function. From the law of addition (|T9|) given in this paper, we learn that 
the function B{z) associated with the unknown G{z) is given by 

B{z) = ^ + zb{z) (69) 

where 

b{z) = Y.Hz) 

i=l 

Thus, we obtain the simple result that the Green's function 
with J2i fi is determined by solving 

^^+biGiz))G{z)=z (71) 

We are however interested in the normalized sum ipT = J2i ^i- De- 
fine the corresponding Green's function as Gt{z) =< j^tr i — - >= 



(70) 

G{z) associated 
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^/KG{\^z). Rescaling, we find that Gt{z) is determined by 



Gr(z) 



]_^..Gt{z) 



Griz) = z (72) 



Note that the argument of bi in (J7%) is ^ . For finite K, (^) is hope- 
lessly complicated. However, as i^T oo we see that it simplifies rather 
naturally to 

^ +a'GT{z) = z (73) 



Griz 
with 



1 ^ 



^'^j^Em (74) 



Solving this quadratic equation, we find immediately that 



Griz) = ^{z- Vz^ - 4a2) (75) 



and thus Wigner's semi-circle law for the density p(yu) = 2^^^'^'^ ~ 

What is 6i(0)? We note that with Gi{z) -* \ + ^ + for large z it 

is easy to show that 6i(0) = af. Thus, not only do we obtain the Wigner 
semi-circle law, we learn that cr^ = -j^Ya cr'f- We should remark that in 
general af is not directly related to the width of the spectrum of eigenvalues 
of ifi, as the reader can easily check by using the explicit formulas given in 
the footnote [27]. 

While we have proved that the density of eigenvalues of (fx satisfies the 
semi-circle law we cannot yet conclude that the probability distribution of ip 
is Gaussian. The reason is that for matrices in the Wigner class, even when 
the distribution is not Gaussian, the corresponding density of eigenvalues 
still satisfies the semi-circle law, as is well known. (For a simple proof based 
on a renormalization group inspired approach, see [12].) For the trace class, 
on the other hand, it is true that if the density satisfies the semi-circle law, 
then the distribution of the matrices is indeed Gaussian [5]. However, we 
see no reason that (fx would belong to the trace class, and not to a more 
involved class of probability distributions such as those discussed in [12]. 

As mentioned above, it is not difficult to extend one of the usual proofs 
of the central limit theorem to the case of matrices. The distribution for the 
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normalized sum matrix (p = is given by (here we omit the subscript 

T) 

(76) 

The integral over can be done in the large K limit: 

(77) 

where we have defined = J dipiP{ipi)j^trip'^. Reexponentiating and in- 
tegrating over t we obtain the desired result that P{(p) is proportional to 



JV i 2 
6 4^2 ^ . 
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Figure Captions 

Fig 1. Feynman rules: (a) quark propagator, (b) gluon propagator, (c) 

quark gluon vertex, (d) gluon interaction, illustrated here with a g(fi^ vertex. 

Fig 2. Quark propagator and one-particle- irreducible self energy. 

Fig 3. Quark self energy: the gluons shown explicitly are all of type 1. There 

are of course also type 2 gluons inside the quark propagator G. 

Fig 4. A class of diagrams not included in (3d). 
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